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ORIGINAL ARTICLE 

A Mechanistic Model of Drug-Induced Liver Injury Aids 
the Interpretation of Elevated Liver Transaminase Levels 
in a Phase I Clinical Trial 

BA Howell 1 , SQ Siler 1 , LKM Shoda 1 ,YYang 1 , JLWoodheacTand PBWatkins 12 

Entolimod (CBLB502) is a Toll-like receptor 5 agonist in development as a single-dose countermeasure against total body 
irradiation. Efficacy can be assessed from animal studies, but the "Animal Rule" does not apply to safety assessment. Marked 
elevations of serum aminotransferases (exceeding 1,000 IU/I) were observed in some human subjects receiving Entolimod in a 
safety study, threatening its continued development.The percentage of total hepatocytes undergoing necrosis in these subjects 
was estimated using a mechanistic, multiscale, mathematical model (DILIsym).The simulations suggested that no subject in 
the safety study experienced more than a modest loss of hepatocytes (<5%), which was comparable to estimates from a study 
of healthy volunteers receiving treatment with heparins. The predicted hepatocyte loss with Entolimod was lower than that 
required to cause liver dysfunction or that is routinely excised from volunteers donating for autologous liver transplantation 
and did not likely represent a serious health risk. 
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INTRODUCTION 

Entolimod (CBLB502) is a Toll-like receptor 5 agonist cur- 
rently in development by Cleveland BioLabs (http://www.cbio- 
labs.com) as a single-dose therapy for exposure to potentially 
lethal whole body irradiation during a nuclear disaster. Dem- 
onstrating efficacy of Entolimod in humans is not feasible, and 
the path to regulatory approval therefore comes under the US 
Food and Drug Administration's (FDA) "Animal rule." 1 The ani- 
mal rule removes the burden of proving treatment efficacy in 
humans. The efficacy of Entolimod, which is a recombinant 
protein injected subcutaneously, has been demonstrated in 
multiple animal species, including mice and nonhuman pri- 
mates. 2 However, the animal rule does not remove the burden 
of proving safety in humans. Accordingly, studies were con- 
ducted in adult healthy volunteers with the sole purpose of 
demonstrating the safety of Entolimod. Although most of the 
subjects tolerated the treatment well in these studies, some 
subjects experienced markedly elevated serum alanine ami- 
notransferase (ALT) and aspartate aminotransferase (AST) 
levels, raising questions about the risk of continued safety 
testing necessary to attain regulatory approval of Entolimod. 

Clinical symptoms were monitored and biomarker measure- 
ments were taken to assess the underlying health status of the 
healthy volunteers. Interpretation of the ALT/AST results was 
not straightforward. Much effort has been put into the devel- 
opment of novel biomarkers aimed at increasing sensitivity, 
specificity, and thus reducing ambiguity. 34 However, many his- 
torical biomarkers of injury or function have remained the gold 
standards for safety assessments in the interim. Liver trans- 
aminase levels in serum are one prime example. ALT and AST 
have, for decades, been the biomarkers of choice for detect- 
ing liver injury (hepatocyte death) in humans. Yet, there are 



relatively little data in humans that link ALT and AST to the loss 
of liver mass. Furthermore, ALT and AST are not entirely liver- 
specific, providing an additional element of uncertainty. As the 
scientific community continues to use historic endpoints such 
as ALT and AST while developing newer ones, the need for 
the interpretation of these safety signals still exists. 

Mechanistic, mathematical models of biological systems 
are ideally suited for biomarker interpretation. They can pro- 
vide the biological context that gives rise to endpoints such 
as ALT and AST. While the application of such models in 
clinical or preclinical situations is still in its nascent state, 
examples of successful application do exist. 5 6 Furthermore, 
the FDA specifically refers to this approach as an important 
component of the modernization of toxicology moving for- 
ward/The chief advantage of using mechanistic models to 
better understand biomarkers is the ability to account for 
all the information inherently included in complicated data 
sets. 

In this report, a mechanistic model of drug-induced liver 
injury (DILI) (DILIsym) was used to aid in the interpretation 
of the liver enzyme elevations observed in healthy volunteers 
receiving Entolimod. This was accomplished through a two- 
phase process. First, DILIsym was optimized to produce the 
liver enzyme profiles observed in healthy volunteers adminis- 
tered Entolimod. Second, the underlying level of hepatocyte 
loss predicted by the model to give such profiles was assessed 
for its relevance to the safety of the healthy volunteers. 

RESULTS 

Baseline human simulations 

Figure 1a shows serum ALT time course data from 100 
healthy volunteers who received a single dose of Entolimod. 
Figure 1b shows three simulated ALT time courses that 
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correspond to the maximum peak ALT observed, the upper 
95th percentile peak ALT observed, and the median peak 
ALT observed (statistical descriptors for measured values 
provided in Supplementary Table S3). The optimization of 
the ALT time courses led to reasonable similarity between 
the simulated ALT time courses and the clinical observations. 
Supplementary Figure S1 shows the AST time courses that 
resulted from the necrosis levels necessary for the ALT data 
fitting. 

In Figure 1c, the peak amount of hepatocyte loss nec- 
essary in the baseline human to generate the maximum 
observed ALT and AST elevations, as well as the 95th percen- 
tile and median elevations, (Figure 1 and Supplementary 
Figure S1) is shown. For most of the healthy volunteers who 
received Entolimod (i.e., the median volunteer), the predicted 
hepatocyte loss was less than 0.5% of their baseline val- 
ues. The upper 95th percentile loss value predicted was 1%. 
The maximum responder with an ALT peak of greater than 
1,000U/I was predicted to have lost a maximum of 3.5% of 
their hepatocytes. 

The time required for the baseline human to regenerate the 
lost hepatocytes was also assessed with DILIsym version 1 A 
(Supplementary Figure S2a). The liver regeneration time 
required to reach 99% of baseline hepatocyte levels for the 
maximum baseline responder (3.5% liver loss) was ~3 weeks. 
The time required for the upper 95th percentile responder to 
recover half the predicted hepatocyte loss (1% liver loss) was 
around 2 weeks. For both cases, the time required for liver 
regeneration was inconsequential, as the level of liver loss 
could probably never be experimentally distinguished from a 
baseline liver. 

Human population sample simulations 

The SimPops constructed from the parameters listed in 
Supplementary Table S1 was composed of 310 simu- 
lated humans. Various levels of necrosis were simulated 
in each individual. An example of this process is shown in 



Supplementary Figure S3. At each level of necrosis, simu- 
lated humans responded with an array of ALT time-course 
profiles. Supplementary Figure S3 includes three exemplar 
induced necrosis levels. Upon review of Supplementary Fig- 
ure S3, it is apparent that simulated variations occur in the 
peak value of ALT but not in the time to reach the peak value. 
This is by design, as the clinical observations from Entolimod 
displayed extremely consistent time-to-peak profiles with dif- 
ferent peak levels. 

A large number of induced necrosis levels were assem- 
bled to produce Figure 2. Figure 2 serves as a look-up table 
(hepatocyte loss vs. peak ALT) for ALT time-course profiles 
that match those exhibited by the Entolimod responders. 
A maximum observed ALT response along the x axis of 
1 ,001 -1 , 1 00 U/l corresponds to the predicted hepatocyte loss 
levels ranging from 2.6 to 4.9% (upper 95th percentile ALT 
response (301-400 U/l) to a range of 0.8-1.7% hepatocyte 
loss, median response to 0-0.5% hepatocyte loss). Supple- 
mentary Table S2 includes the ranges associated with each 
peak ALT category shown on the x axis of Figure 2. 

The SimPops was also used to assess variability in liver 
regeneration time after hepatocyte loss. The time required 
for liver regeneration in the simulated individuals with ALT 
levels corresponding to the maximum observed value in the 
study ranged from 2 to 9 weeks with a median value of 3 
weeks. The results are displayed in Supplementary Figure 
S2b as a shaded region, which encompasses all simulated 
observations. 

DISCUSSION 

Peak serum levels of liver enzymes observed in a few sub- 
jects taking Entolimod raised safety concerns, but the extent 
of liver injury that had occurred, and hence the safety risk to 
these subjects, was unclear. The levels of predicted hepa- 
tocyte loss shown in Figures 1c and 2 can generally be 
described as "low" for individuals receiving the single dose 
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Figure 1 Simulation results for baseline human compared to observed responses in a clinical trial of 100 normal healthy volunteers receiving 
Entolimod. (a) Alanine aminotransferase (ALT) profiles for healthy volunteers are shown with the blue continuous lines, (b) The maximum 
ALT value observed (black line) was used along with the ALT time-course dynamics to optimize DILIsym model simulations to the maximum 
observation (red line). The same procedure was followed for the upper 95th percentile ALT observation and the median ALT observation. 
Simulation results are shown in red. ALT time courses are shown in black, (c) The maximum level of hepatocyte loss after Entolimod dosing 
relative to the number of healthy hepatocytes in adults was predicted with DILIsym using simulations corresponding to the maximum, 95th 
percentile, and median peak ALT responses measured in the clinical trial. 
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Figure 2 Range of predicted hepatocyte loss (percentage of 
baseline) as a function of observed peak alanine aminotransferase 
(ALT) range in healthy volunteers. The range of hepatocyte loss 
was predicted using version 1 A and SimPops (in silico population 
samples) specifically designed to replicate physiologically relevant 
variability in ALT response from necrosis of hepatocytes. The dashed 
lines indicate the predicted range of hepatocyte loss associated with 
a peak ALT range consistent with the maximum ALT observed in 
normal healthy volunteers dosed with Entolimod. 

of Entolimod. Even the estimated range of hepatocyte loss 
in the maximum ALT responder (2.6-4.9%) intuitively seems 
modest given the enormous functional reserve of the liver 
and its remarkable regenerative capabilities. To put this esti- 
mated hepatocyte loss into context, three types of published 
studies were used or reviewed below. 

Comparing Entolimod responses to heparin responses 

Entolimod responses were put into context regarding safety 
by comparing them to the responses from another agent 
that is generally regarded as safe for the liver. Harrill et al. 
analyzed serial serum samples obtained from healthy vol- 
unteers receiving several types of heparins. 3 ALT elevations 
were observed, and experimental serum biomarkers (rises 
in high-mobility group box 1 and glutamine dehydrogenase 
with no rise in caspase-cleaved fragments of Keratin-18) sup- 
ported necrosis as the underlying hepatocyte death process. 
Data from the two volunteers who experienced the highest 
observed ALT elevations during heparin treatment are shown 
in Figure 3a. ALT dynamics from heparins are very different 
when compared to Entolimod responses (Figure 1a). Figure 
3a also shows two solid lines indicating DILIsym simulation 
results. The same parameter values used for Entolimod simu- 
lations (Table 1) were used to produce the ALT simulations 
that match the time-course data from the heparin study. The 
primary difference was the time for onset of injury, which was 
markedly lengthened for the heparin simulations. Figure 3b 
shows that baseline predictions of liver loss in the two highest 
ALT responders in the heparin study were 4.4% and 2.7%. 
The level of liver loss predicted for the heparin responders 
was on the same order as the maximum values predicted for 
the Entolimod responders. Since heparins are regarded as 
safe for the liver, 38 this comparison supports the conclusion 
that the hepatocyte loss predicted with Entolimod can occur 
with other therapeutic agents without posing a significant 
safety risk to those treated. 
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Figure 3 Clinical data and simulation results related to alanine 
aminotransferase (ALT) elevations observed in healthy volunteers 
receiving various heparins. 3 (a) The two volunteers who exhibited 
the highest observed ALT elevations (data points) in a previously 
published clinical study utilizing heparins, 3 and DILIsym simulation 
results (solid lines) from a model optimization designed to fit the 
heparin observations, (b) The predicted maximum percentage of 
hepatocytes lost for the two heparin responders using slow injury 
onset and fast ALT transport (see Figure 5). Data reprinted with 
permission. 

Comparing predicted hepatocyte loss levels to historical 
data 

A direct comparison between predicted levels of liver loss for 
Entolimod and levels of liver loss that lead to clinically sig- 
nificant symptoms of liver dysfunction was conducted using 
published data. 9 Figure 4 shows prothrombin time and bili- 
rubin as a function of hepatocyte volume fraction for adults. 
Also, illustrated in Figure 4 are the levels of hepatocyte vol- 
ume fraction predicted for the Entolimod responders (note 
that not all liver volume is composed of hepatocytes, so the 
levels indicated on the x axis are lower than the predicted 
percentage of hepatocytes lost). The levels of loss implicated 
for Entolimod are below the threshold for detectable liver dys- 
function (Figure 4a,b), which is consistent with expectation, 
since there were no elevations in serum protime observed 
in the clinical study. Transient and minor bilirubin elevations 
were observed in some of the volunteers, but there was no 



www.nature.com/psp 



DILIsym Aids in Liver Safety Assessment in a Phase I Clinical Trial 

Howell eta/. 



Table 1 Key DILIsym version 1A parameters related to the Entolimod application 



DILIsym parameter (units) 



DILIsym parameter 
description 



Parameter 
value used 



Lower 
bound 



Upper 
bound 



HC_necr_ALT_release (units/1 e9 hepatocytes) Hepatocellular ALT content 

Liver_ALT_h, interm_ALT_h, serum_ALT_h b (hour) ALT half-life 

lnterm_ALT_serum_ALT_rate, liver_ALT_interm_ALT_rate c (1/hour) ALT transport rate constant 

HC_necr_AST_release (units/1 e9 hepatocytes) Hepatocellular AST content 

Liver_AST_h, interm_AST_h, serum_AST_h b (hour) AST half-life 

lnterm_AST_serum_AST_rate, liver_AST_interm_AST_rate c (1/hour) AST transport rate constant 



530 a (ref. 5) 250 (ref. 22) 77,242 (ref. 23) 



47 (ref. 24) 37 

2 d 0.08 (ref. 17) 

1,079 (ref. 25) 424 (ref. 22) 

18 (ref. 24) 12 

2 d 0.08 (ref. 25) 



Pro_regen_coeff (dimensionless) 



Hepatocyte proliferation rate constant 7 (refs. 9,26,27) 



57 
2.4 
1,358 
24 
2.4 
10 



ALT, alanine aminotransferase; AST, aspartate aminotransferase. 

a Vastly different values for this parameter are reported, but Remien etal. (201 2) 5 reports a value consistent with the AST content reported in Rej etal. (1978) 25 
and the ALT content reported in Boyd (1983), 22 so this value was chosen as the default value for this application. b Although the model gives users the option 
of adjusting the half-lives for each ALT/AST compartment individually, the values were kept equal for this application. This is partly due to the knowledge gap 
concerning the primary clearance mechanism for ALT and AST from the plasma. c Two parameters are necessary to describe the transport between the three 
ALT/AST compartments, but the values are generally identical for those parameters within DILIsym. Including the two parameters separately gives the user the 
option of altering this setup. d The value of the parameter used for this application was based on the clinical data provided for Entolimod, where the rapid rise in 
ALT/AST observed specified that the value be large compared to the baseline human in DILIsym version 1 A to recapitulate the data. This value was also used 
for the heparin optimization shown in Figure 3. The lower bound value of 0.08 leads to ALT peak times consistent with clinical acetaminophen data when the 
ALT half-lives are specified as listed here. 



correlation between peak serum bilirubin levels and peak 
ALT values (Supplementary Figure S4), suggesting a 
mechanism independent of liver cell death. The subject with 
the highest ALT elevation had a serum bilirubin peak that 
never reached or exceeded 1 mg/dl, which is below the upper 
limit of normal (ULN). Furthermore, bilirubin elevations never 
exceeded 2 X ULN in any volunteer, the typical cut off for clin- 
ical concern. 10 Accordingly, these observations support the 
conclusion that the maximum hepatocyte loss experienced 
in the clinical study was far short of the extent necessary to 
detectably interfere with liver function. 

Liver excision in healthy adults 

One practical way to provide context to liver loss predictions is 
through comparisons to liver excision levels in healthy donors 
for the purpose of autologous transplantation. For healthy 
adult donors participating in adult-to-pediatric liver transplan- 
tation, removal of the left lateral segment of the liver, which 
constitutes roughly 20% of the total liver volume, is consid- 
ered of "no concern with regard to the donor's remnant size." 11 
Furthermore, when left or right lobes are used for adult-to- 
adult living donor liver transplantation, roughly 40-60% of the 
living donor's liver volume is removed, respectively. 12 A clean 
excision of liver tissue for transplantation is different than DILI 
in many ways. Nonetheless, the large masses of liver loss 
occurring during transplantation provide some context for 
the loss of 5% predicted for the most extreme responder to 
Entolimod in the healthy volunteer study. It is also interesting 
to note that such donations of liver by healthy volunteers is 
considered ethical, which supports the idea that continued 
evaluation of Entolimod in healthy volunteers may also be 
ethical with appropriate informed consent. 

Interpreting liver enzyme signals in the context of safety 

In this example, the liver enzyme elevations from healthy vol- 
unteers (Figure 1 and Supplementary Figure S1) appeared 
rapidly and disappeared rapidly. However, if another com- 
pound gave a similar peak ALT but a different ALT kinetic 
profile, it could alter the potential conclusions from tools such 
as DILIsym. To better illustrate the factors that should be 



considered when examining liver enzyme kinetics, DILIsym 
was used to produce four simulated ALT profiles (Figure 5a). 
Two rates or parameter values were set for each of the four 
scenarios shown, and necrosis was then induced ("tuned") 
until the peak ALT value reached 1 ,000 U/l. The two rates set 
were (i) rate of injury onset, which is dictated by the time 
course for drug exposure and injury progression and (ii) the 
rate of ALT transport from necrosing cells into plasma, which 
could be affected by the mode of cell death. Note the clear 
differences in kinetics between the four curves despite simi- 
lar peak levels. Figure 5b contains the four peak hepatocyte 
loss predictions corresponding to the four ALT time courses 
in Figure 5a. Each ALT time course demonstrates a unique 
profile and a unique level of liver loss. Accordingly, some key 
factors to consider when examining liver enzyme elevations 
can be ascertained from the figure: 

1 . ALT peak height this is an important component but 
clearly not sufficient to distinguish between levels of 
liver loss, as Figure 5b demonstrates. 

2. Area under the curve: Figure 5a shows that both curves 
for fast onset have lower area under the curve (AUC) 
values than the curves for slow onset. While AUC can 
thus be correlated to hepatocyte loss, the predicted loss 
levels are also different between the fast curves, and 
between the slow curves, despite similar AUC values. 
Therefore, AUC is an important component for the inter- 
pretation of the liver enzyme time courses but is not suf- 
ficient alone. 

3. Time between the initiation of liver enzyme elevations 
and the point at which the maximum observed value 
is reached: This time lapse depends on both the onset 
and transport rates and is the final component neces- 
sary to estimate hepatocyte loss from the liver enzyme 
profiles. Note how the injury extent predicted in Figure 
5b gradually increases as this time lapse increases in 
Figure 5a. 

These observations and Figure 5a, b demonstrate the 
value of considering the full depth of information obtained 
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Figure 4 Symptoms of hepatocyte loss relative to the magnitude 
of loss, (a) Relationship between the maximum prolongation of 
prothrombin time in the first 10 days after acetaminophen overdose 
and the hepatocyte volume fraction in 76 patients. 9 (b) Relationship 
between the maximum plasma bilirubin in the first 10 days after 
acetaminophen overdose and the hepatocyte volume fraction in 75 
patients. 9 The dashed line represents a significant clinical bilirubin 
elevation of 3mg/dl. Figures reprinted with permission. The range of 
predicted hepatocyte volume fraction for the volunteers exposed to 
Entolimod is also indicated in both figures. 



from time-course profiles of biomarkers. They also highlight 
the importance of the exact time-course profiles observed 
for Entolimod. Had the Entolimod healthy volunteers exhib- 
ited liver enzyme profiles similar to the ones seen for hepa- 
rin responders with maximum observations near 1,00011/1, 
the liver loss predictions would have been roughly threefold 
higher than the current predictions, supporting a greater 
safety concern. 
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Figure 5 Alternate parameter combinations yielding simulated 
alanine aminotransferase (ALT) profiles with similar peak values but 
different time-course dynamics and the simulated peak hepatocyte 
loss associated with each ALT profile, (a) Four simulated ALT 
profiles are shown. The ALT profile most similar to the profiles 
from Entolimod responders (solid black line) was generated with 
fast injury onset and fast transport of ALT out of necrotic cells into 
plasma. Fast injury onset and slow transport into the plasma (dotted 
black line) is also shown. The profile most similar to the heparin 
observations shown in Figure 5 was simulated with slow injury 
onset and fast ALT transport (solid gray line). An alternative profile 
that is also consistent with the heparin observations includes slow 
injury onset and slow ALT transport (dashed gray line), (b) Simulated 
peak percentages of hepatocytes lost from the DILIsym model 
corresponding to the profiles shown in a. 

While the connection between hepatocyte death and liver 
enzyme leakage, as well as the use of time-course profile 
information to interpret liver enzyme elevations may seem 
intuitive in the context of this example, it must be noted that 
the current FDA guidelines surrounding the assessment of 
DILI in clinical trials lack any mechanistic or physiological 
basis. 10 For example, the FDA guidance recommendation to 
stop investigational drug treatment if a subject experiences a 
rise in ALT exceeding 8 X ULN on one occasion, or exceed- 
ing 5 X ULN for more than 2 weeks 10 is empirical. The meth- 
ods for analyzing liver enzyme elevations outlined herein, 
together with estimations of the relative contribution of necro- 
sis and apoptosis that can now be obtained from mechanis- 
tic biomarkers, could represent a significant advancement in 
regulatory science in the future. 
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Limitations 

Mathematical models include approximations that limit their 
accuracy. Clinical experience should weigh heavily on con- 
clusions regarding safety. Furthermore, the selection of 
model parameters to create population samples for variability 
predictions could lead to parameter value combinations that 
are not biologically reasonable. In this example, a small num- 
ber (four) of seemingly independent parameters were used to 
create the SimPops, and those parameters were directly tied 
to observed measurements. 

In this case study, the goal was to answer one pri- 
mary question: how much liver loss would have to occur 
to result in the given ALT/AST profiles? As a result, the 
specific mechanisms of Entolimod induced enzyme eleva- 
tions were neglected. Specific pathways for inducing cell 
death by Entolimod could also be explored with DILIsym 
if the appropriate time and data were available. Necrosis 
was assumed to have occurred as the primary mode of cell 
death. While this assumption does seemingly provide the 
most conservative approach, alternative hypotheses, such 
as apoptosis, could lead to different conclusions regard- 
ing hepatocyte loss. Had serum been available from the 
Entolimod studies, it would have been helpful to assay 
serum biomarkers capable of distinguishing cell death by 
necrosis and apoptosis. 3 In the future, application of these 
biomarkers in clinical trials would provide very valuable 
input for modeling exercises aimed at improving liver safety 
assessment. 

Finally, this case study was a retrospective analysis of 
data collected during two healthy volunteer studies involv- 
ing Entolimod. The simulations and associated conclusions 
were meant to aid in determining the implications of those 
liver signals for the healthy volunteers involved but were not 
meant to prospectively predict or be applied to the safety 
of different healthy volunteers or patients chosen for future 
studies. This is especially true in light of the lack of mechanis- 
tic detail included for the simulation of liver enzyme increases 
in plasma by Entolimod. 

Summary and conclusions 

In conclusion, DILIsym version 1A predicted hepatocyte 
loss levels of less than 1% for most of the healthy vol- 
unteers receiving Entolimod in the completed clinical tri- 
als. The maximum responder's loss was predicted to fall 
between 2.6% and 4.9%. In light of heparin data, historical 
patient biopsy information, and liver excision practices, the 
amount of hepatocyte loss experienced by the volunteers 
was unlikely to have represented a significant health threat. 
DILIsym version 1A also demonstrated the importance of 
considering the dynamics of liver enzyme data. The three 
most important factors to consider include the maximum 
value observed, the AUC, and the time elapsing between 
the start of enzyme elevation and when the maximum value 
was observed. 

METHODS 

Safety studies in healthy volunteers 

The sponsor of Entolimod conducted two safety studies in 
healthy adult volunteers. They included 50 individuals (dosing 



range of less than 5 ug up to 50 ug as a single dose) and 
100 individuals (dosing range of 25-35 ug). The results pre- 
sented in Figure 1 are constrained to the larger study, where 
the most significant ALT elevations were observed (patient 
demographics in Supplementary Table S4. 

The primary endpoints analyzed were biomarkers (ALT, 
AST, and serum bilirubin) measured from serial serum sam- 
ples. The focus of the modeling effort was ALT because liver 
safety concerns were based on these measurements. The 
main goal of this exercise was to use the aminotransferase 
data to estimate the percentage of hepatocytes lost in the 
livers of the healthy volunteers following the dose of Ento- 
limod. The mode of induced hepatocyte loss by Entolimod, 
which may be related to its properties as a Toll-like receptor 
5 agonist, were not included in the simulations, nor were the 
pharmacokinetics of the drug represented. Rather, this study 
focused on the modeled relationships between hepatocyte 
necrosis, the resulting release of hepatocyte ALT and AST, 
and its distribution to and clearance from the plasma com- 
partment. With this representation, simulations were used to 
predict the level of hepatocyte death that would lead to the 
observed elevations in ALT and AST. 

DILIsym version 1A 

The DILI-sim Initiative is a collaboration between the Ham- 
ner Institutes for Health Sciences and the pharmaceutical 
industry. The goals of the DILI-sim Initiative include develop- 
ing DILIsym software, a predictive, mechanistic, mathemati- 
cal model of DILI (http://www.DILIsym.com). The method for 
model design is best described as "middle out." The "middle 
out" approach involves starting at the organ level (liver in this 
case) and working down to the molecular level or up to the 
organism level when necessary. 13 The model is organized 
into various smaller submodels, but all submodels are math- 
ematically integrated to simulate an organism level response. 
Supplementary Figure S5 shows several of the key compo- 
nents of DILIsym version 1 A and how they exist as submod- 
els but also interact. More information on DILIsym version 1 A 
is available through several previously published articles. 6 ' 1415 
DILIsym is available to industry through membership in the 
DILI-sim Initiative. Academic and nonprofit groups may 
access the model by contacting the Hamner Institutes for 
Health Sciences or the corresponding author. 

ALT and AST submodels of DILIsym version 1 A 

This exercise used the hepatocellular life-cycle and bio- 
marker components of DILIsym version 1 A (Supplementary 
Figure S5, more detail including equations in Supplemen- 
tary Materials and Methods). The ALT submodel was the 
principle focus. The AST submodel, which contains an identi- 
cal structure with different parameter values, was also used. 
The ALT submodel is shown in Figure 6. The necrosis of 
hepatocytes leads to ALT release. 16 The mass of ALT per 
cell is specified. Three compartments (modeled with ordinary 
differential equations) make up the ALT model. In Figure 6, 
the compartments are labeled "Hepatocyte," "Intermediate," 
and "Circulating." The "Hepatocyte" compartment repre- 
sents temporary association of ALT with cells as necrosis 
occurs. The "Circulating" compartment represents systemic 
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Figure 6 Diagram showing the structure of the alanine transaminase 
(ALT) biomarker submodel within DILIsym version 1 A. 



circulation. The "Intermediate" compartment was required 
to optimize observed serum ALT profiles with other model 
hepatotoxicants and may represent the ALT in interstitial 
space, the Space of Disse, or possibly the liver sinusoids. 
As shown in Figure 6, clearance of ALT occurs in all three 
compartments. This structure was implemented based on 
ALT dynamics observed in response to several hepatotoxic 
compounds. 417 " 19 Note that the clearance of ALT occurs in all 
three compartments, as the exact mechanism and location 
of liver enzyme clearance is not yet elucidated. 2021 

Table 1 shows several of the key parameter values for the 
ALT and AST submodels used to match the ALT responses 
observed in healthy volunteers given the single dose of Ento- 
limod. 5 9 17 22-27 The values used for the response in the base- 
line (average) human, as well as the range of values derived 
from the literature, are shown. The time to reach peak ALT/ 
AST values was much faster with Entolimod than with previ- 
ous exemplar hepatotoxic compounds modeled (e.g., acet- 
aminophen, furosemide, aflatoxin). 6 ' 1428 This is most likely 
due to the different mechanisms of toxicity for the previous 
exemplars compared to Entolimod. Based on the ALT profiles 
and assuming the ALT and AST were hepatocyte-derived, 16 
hepatocyte death and liver enzyme release appeared rapid. 
As a result, the transport dynamics were adjusted to match 
the observed profiles. 

Baseline human simulations: ALT and AST 

To assess the level of liver injury necessary to produce the 
liver enzyme responses seen with Entolimod, hepatocyte 
necrosis was simulated in DILIsym. Hepatocyte necrosis was 
chosen as the mode of cell death because the sponsor of 
Entolimod lacked any specific biomarker data implicating a 
specific mode of cell death in the human studies (such as 
serum high-mobility group box 1 assessments). 34 The rapid 
nature of the liver enzyme responses observed in serum 
also made necrosis seem most plausible. The health con- 
sequences of necrosis, which could include activation of an 
innate immune response, are potentially more dire compared 
to hepatocyte death by apoptosis. Necrosis was thus the logi- 
cal and conservative choice for this analysis. 



The mechanisms operative before the onset of necro- 
sis (e.g., pharmacokinetics, adenosine triphosphate (ATP) 
depletion, oxidative stress, etc.) were not considered herein. 
The method for simulating necrosis was a direct necrosis 
signal. The amount of necrosis was the starting place for 
the optimization in this example. Necrosis was tuned to 
achieve similar time-course profiles as observed for the 
maximum, upper 95th percentile, and median peak ALT 
level responses to Entolimod. The overall maximum, upper 
95th percentile, and median response levels were calcu- 
lated from the peak ALT values from the entire group of trial 
participants. The parameter values shown in Table 1 were 
used and the amount of necrosis was adjusted until simu- 
lated ALT profiles were consistent with the clinical data. The 
magnitude of the ALT responses were the optimization data 
set. Simultaneously, AST levels were simulated with DILI- 
sym to ensure a reasonable prediction of the magnitudes 
of both enzymes. The AST responses were the confirma- 
tion data set. Next, the simulated percentage of hepatocytes 
lost corresponding to each of the three enzyme profiles was 
inferred from the simulation results. 

Assessing variability with a human population sample 

The simulation process above assumes that the interpa- 
tient variation in serum ALT profile only reflects variation 
in the extent of necrosis. In order to account for potential 
intersubject differences in the model parameters, a popu- 
lation sample (SimPops) was constructed using param- 
eters that would cause variability in ALT response from a 
given level of necrosis. A parameter related to liver regen- 
eration capacity was also used, since this could theoreti- 
cally affect the peak loss of hepatocytes while not affecting 
the serum ALT profile. A population sample with 310 
simulated humans was created and used in the final analysis. 
Various levels of necrosis were then induced for each simu- 
lated human and the simulated variability in ALT response was 
measured. This information was used to construct Figure 2. 

Details related to the construction of the SimPops, includ- 
ing the specific parameters varied for this application, can be 
found in the Supplementary Materials and Methods. Fur- 
ther details on the methods are also in the supporting mate- 
rials of Howell etal. (201 2) 14 and Woodhead etal. (201 2). 6 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

S Healthy volunteer studies were conducted to 
assess the safety of Entolimod, a radiation sick- 
ness countermeasure in development. Interpre- 
tation of the results with regard to liver safety 
were not straightforward. 

WHAT QUESTIONS DID THIS STUDY ADDRESS? 

S How much hepatocyte necrosis likely occurred 
in the healthy volunteers receiving Entolimod to 
result in the measured ALT/AST profiles? 



WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

y The amount of hepatocyte loss experienced by 
the volunteers receiving Entolimod was unlikely 
to have represented a significant health threat. 
The importance of considering the dynamics of 
liver enzyme data was also highlighted. 

HOW MIGHT THIS CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

S Mechanistic, mathematical models of biological 
systems are ideally suited for biomarker inter- 
pretation and can provide the biological context 
that gives rise to endpoints such as ALT and 
AST. Their continued use and adoption could 
prove valuable to drug developers, physicians, 
and regulators when concerns over safety are 
raised. 
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